Density profiles and surface tensions of polymers near colloidal surfaces 



A. A. Louis\ P.G. Bolhuis^, E.J. Meijer^ and J. P. Hansen^ 
^Department of Chemisty, Lensfield Rd, Cambridge CB2 lEW, UK 
^Department of Chemical Engineering, University of Amsterdam, Nieuwe Achtergracht 166, NL-1018 WV Amsterdam, 

Netherlands. 
(February 1, 2008) 

The surface tension of interacting polymers in a good solvent is calculated theoretically and by 
computer simulations for a planar wall geometry and for the insertion of a single colloidal hard- 
sphere. This is achieved for the planar wall and for the larger spheres by an adsorption method, and 
for smaller spheres by a direct insertion technique. Results for the dilute and semi-dilute regimes 
are compared to results for ideal polymers, the Asakura-Oosawa penetrable-sphere model, and to 
integral equations, scaling and renormalization group theories. The largest relative changes with 
density are found in the dilute regime, so that theories based on non-interacting polymers rapidly 
break down. A recently developed "soft colloid" approach to polymer-colloid mixtures is shown to 
correctly describe the one-body insertion free-energy and the related surface tension. 



I. INTRODUCTION 

Binary mixtures of polymers and colloidal particles in 
various solvents are the focus of sustained experimen- 
tal and theoretical efforts, both because they constitute 
a challenging problem in Statistical Mechanics of "soft 
matter" , and because of their technological importance 
in many industrial processes. One of the most striking as- 
pects of polymer-colloid mixtures, namely the depletion 
interaction between colloids induced by non-adsorbing 
polymer was recognized nearly 50 years ago 0|. More 
recently, the importance of the polymer depletant in de- 
termining the phase behavior of the mixtures was realized 
1^, and much recent experimental work was devoted to 
the phase diagram (3[j^, structure interfaces d), 

and to the direct measurement of the effective interac- 
tions |p|-[rT|]. On the theoretical side, most efforts have 
concentrated on impenetrable spherical colloids, while 
various models and theoretical techniques have been in- 
vestigated for the description of the non-adsorbing poly- 
mer coils. The models include non-interacting (ideal) 
polymers | |l | ,|l2| , p^ , polymers represented as penetrable- 
spheres [p~4pp7|, and interacting polymers coarse-grained 
to the level of "soft colloids" p8| - 2l[ | . Monomer level rep- 
resentations of polymer chains, like the self-avoiding walk 
(SAW) model, appropriate in good solvent, have been 
considered within polymer scaling approaches [p2|-p4[, 
renormalization group (RG) theory [p5| - |29f| , and fluid in- 
tegral equations [ ^J3l| . 

While many effects for the simplest case of colloids 
mixed with non-interacting polymer are quantitatively 
understood, the behavior of the experimentally more 
relevant case of polymers with excluded volume inter- 
actions is at best understood on a qualitative basis; a 
quantitatively reliable theory is still lacking. Clearly, 
to construct such a theory for finite concentrations of 
colloidal particles, one must first understand how in- 
teracting polymer coils distribute themselves around a 



single spherical colloid of radius Re- This problem is 
addressed in the present paper using a combination of 
Monte Carlo (MC) simulations and scaling theories to 
determine the key quantities, which are the monomer or 
center-of-mass (CM) density profiles p(r) of SAW poly- 
mers around a single impenetrable sphere, as well as the 
resulting surface tension. If Rg denotes the radius of gy- 
ration of the polymers, these quantities clearly depend 
on the ratio q = Rg/Rc, which controls the curvature 
effects. The limit q —^ 0, corresponding to a polymer 
solution near a planar wall, will be examined first, be- 
fore considering the case of finite q. The complete the- 
ory for the opposite limit, g 3> 1, will be the subject 
of a future publication, although we show some prelimi- 
nary results here. Throughout this work we focus on the 
dilute and semi-dilute regimes |^,^ of the polymers, 
where the monomer density c is low enough, for detailed 
monomer- monomer correlations to be unimportant; the 
melt regime, where c becomes appreciable, will not be 
treated here. 

The surface tension of a polymer solution surround- 
ing a sphere is macroscopically defined by considering 
the immersion of a single hard colloidal particle into a 
bath of non-adsorbing polymer. Because this immer- 
sion reduces the number of configurations available to the 
polymers, resulting in an entropically induced depletion 
layer around the colloid, there is a free energy cost Fi 
for adding a single sphere to the polymer solution which 
naturally splits into volume and surface terms: 



i^i =n(p)-7ri?3+4^i?2^,(p). 



(1) 



The first term in Eq. (|l|), describes the reversible work 
needed to create a cavity of radius Rc in the polymer 
solution. Since the osmotic pressure Ii{p) of a polymer 
solution in the dilute and semi-dilute regimes is quanti- 
tatively known as a function of polymer concentration p 



1 



from RG calculations this volume term is well un- 
derstood. The problem of a quantitative description of a 
single colloid in a polymer solution thus reduces to un- 
derstanding the second term, which defines the surface 
tension "fsip), i-c the free-energy per unit area that is 
directly related to the creation of the depletion layer. It 
is customary to relate the surface tension 7s (p) around 
a sphere to the surface tension 7^ (p) near a planar wall, 
by expanding in powers of the ratio q =^ Rg/Rc'. 



(2) 



which is expected to be most useful when q is not too 
large. The coefficients control the curvature cor- 

rections. They are analogous to the Tolman corrections 
in the macroscopic case [^jsi . 



The paper is organized as follows: The case of a single 
plate or hard wall immersed in a polymer solution is dis- 
cussed in section II, where we report results for density 
profiles p{z) at various polymer concentrations. These 
density profiles define the reduced adsorption T{p), from 
which the surface tension "fw{p) may ultimately be ex- 
tracted. These considerations are extended to spherical 
colloids in section III, where simulation results for the 
density profiles are reported for size ratios q — 0.67, 1.05 
and 1.68. These data are then used to compute 7s (p) and 
the Ki(p); limiting forms are extracted for the p — > and 
the semi-dilute regimes. The results are compared to the 
theoretical predictions for ideal polymers, for the pene- 
trable sphere model, and wherever applicable, to RG and 
integral equation predictions. The limit of large q, where 
the expansion (^) becomes less useful is also discussed. 
For this limit we also report on some preliminary direct 
simulation results for Fi based on the Widom insertion 
technique js^. Finally we show that the "soft colloid" 
paradigm has the correct thermodynamics of the single 
colloid problem automatically built in. 



II. DENSITY PROFILES AND SURFACE 
TENSION NEAR A SINGLE WALL 

A single hard wall in a bath of non-adsorbing poly- 
mers creates an entropically induced depletion layer be- 
cause the polymers have fewer possible configurations 
near the wall. To calculate these density profiles we 
performed Monte-Carlo simulations of the popular self 
avoiding walk (SAW) model on a cubic lattice. Even 
though this model ignores all chemical details of a real 
polymer system except the excluded volume and poly- 
mer connectivity, it reproduces the scaling behavior and 
many other physical properties of athermal polymer so- 
lutions [^,^ . For N polymers of length L on a lattice 
of M sites, the polymer density is given by p = N/M, 
while the monomer density is c = LN/M. The poly- 
mers are characterized by the radius of gyration, which 
scales as Rg ^ , where v « 0.588 is the Flory exponent 
p2|,p5|j32|. For densities p less than the overlap density 



p* = l/(|7r_Rg) the system is in the dilute regime, while 
for p > p*, and c ^ 1, the system is in the semi-dilute 
regime. We use L = 500 SAW polymers in our simula- 
tions, which are expected to exhibit properties close to 
those corresponding to the scaling limit L ^ 00. Further 
details of the simulation method and the model can be 
found in [p^-pl|. Note that a small correction to these 
results must be applied ]3^ ]! Since our models are all 
athermal, we set /? = l/{kBT) — 1. 

Examples of the depletion layer density profiles near a 
hard wall are depicted in Fig. |l| for a polymer center of 
mass (CM) representation, as well as for a monomer rep- 
resentation. Both profiles have, by definition, the same 
reduced adsorption, defined as [p4|: 



r(p) 



1 din"'' /A) 
P dp. 



h{z)dz, 



(3) 



where D,'^^ /A is the surface excess grand potential per 
unit area A. h{z) = p{z)/ p — 1, with p{z) the CM den- 
sity profile of the polymer coils a distance z from the 
surface. In the monomer representation one should re- 
place p by c = LN/M and h{z) by the monomer profile; 
the two reduced adsorptions are equal and measure the 
reduction in the number of chains near the surface. 

In the low-density limit an RG calculation based on 
a first order e-expansion gives f (0) ~ — 1.074i?g [p7| , ^ , 
which is slightly less than 



f*'' = 2R, 



-1.128i?, 



■9' 



(4) 



the density independent result for an ideal polymer with 
the same size Rg ||2^ (but larger L due to the different 
scaling of the radius of gyration) . This reflects the fact 
that for a given the polymer-polymer interactions 
reduce the size of the depletion layer, an effect which be- 
comes more pronounced with increasing density, see e.g. 
Fig. % 

For the semi-dilute regime, de Gennes has proposed 
an approximate expression for the monomer profile near 
a wall, h„i{z) = tanh^(z/^(p)) — 1, where ^(p) is the 
correlation length or blob size |^^. If we identify ^(p) 
with — r(p) then, as shown in Fig this form provides 
a fairly accurate fit to our simulation results. Since 



C(p) ~ p— /(3-1) 



-0.770 



in the semi-dilute regime. 



this implies that the density profiles should become more 
narrow with increasing density, a trend clearly seen in 
Fig.0. 

From the density profiles of Fig. |l|, we can derive the 
adsorption at different densities using Eq. (||). These 
are shown in Fig. H, together with a simple fit con- 
strained to give the expected p = value, and the cor- 
rect scaling behavior in the semi-dilute regime where 
-f(p) « ^(p) - p-"-'^™, namely 



r(p) = 



-1.074i?„ ( 1 + 7.63— + UM(—f 



-(0.2565) 



(5) 
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Throughout this paper the value of the radius of gyra- 
tion is conventionally chosen as that appropriate for an 
isolated polymer, i.e. Rg — Rg{p = 0). However, as 
the polymer concentration increases, the measured Rg{p) 
will decrease with density as shown in Fig. ^. In the 
semi-dilute regime this scales as Rg{p)/ Rg{p = 0) ^ 
^(2iy~i)/(6i'-2) ^ ^-0.115 HIj^ which decreases much more 
slowly with density than the correlation length ^{p) or 
the relative adsorption T{p). In fact at p/p* = 1, the 
crossover from the dilute to the semi-dilute regimes, T{p) 
has dropped to 59% of its p — > value, while Rg{p) has 
only changed by a few %. The largest rate of relative 
change in the adsorption is therefore found in the di- 
lute regime, suggesting that theories based on the p ^ 
limit may start to break down well before the semi-dilute 
regime is reached. The border between the dilute and 
semi-dilute regimes is not sharp. For the semi-dilute 
regime, the asymptotic forms derived by scaling theories 
appear to be reached at a lower density for T(p) than for 

Rgip) in. 

One route to calculate the surface tensions from the 
density profiles is to use an extension of the Gibbs ad- 
sorption equation to express the surface tension near a 
single wall in terms of the relative adsorption and the 
equation of state: 



dA 



" fdn{p') 



V dp' 



tip')dp'. 



(6) 



The derivation of this equation can be found, for exam- 
ple, in p8| , p4[ . By performing one integration by parts 
w.r.t. density, Eq. (|6|) can also be expressed as: 



7^(p) = -n(p)f(p)+ / n(p') 



dp' 



dp', (7) 



The first term in this equation takes the form of a pres- 
sure times a length. For ideal polymers, where T{p) is 
independent of density [|l],^, this term completely de- 
scribes the surface tension of the depletion layer. It is just 
the (entropic) free energy cost per unit area of creating 
a cavity of volume r{p)A. The second (positive) term 
is therefore only relevant if there are polymer-polymer 
interactions. 

We have previously calculated the equation of state 
for L = 500 and L = 2000 SAW polymers |l9|,0, 
both of which are well described by analytic RG ex- 
pressions Using this for Il{p) together with the 
fit to the relative adsorption from Eq. (H), we can now 
use Eq. (0) to calculate the surface tension of a solu- 
tion of polymers in good solvent near a single wall. Our 
results are shown in Fig ^. In the low density limit 
the surface tension reduces to the same functional form 
as for ideal polymers, i.e. jw{p) = — n(p)r(/9), so that 
limp/p.^o 7u,(p) « 1.074i?gp. Note that for all but the 
lowest densities, the surface tension is considerably larger 
than the ideal polymer result 7^(/3) — 2pRg/^jT:. The 
surface tension for interacting polymers increases more 



rapidly with increasing density both because n(p) in- 
creases faster than T{p) decreases in the first term of 
Eq. (0), and because the second term, which is absent 
for non- interacting polymers, increases with density as 
well. 

Further simplifications occur in the semi-dilute regime. 
For example, when the scaling forms for the osmotic pres- 
sure, n ^ p3i//(3i/-i)^ g^-^^ ^Yic reduced adsorption, 
f ~ p-i^/CS"-!)^ are used in Eq. (Q), then the integral in 
the second term can be easily performed and turns out 
to be exactly half the first term, a result that is indepen- 
dent of the value of the exponent v. The surface tension 
therefore takes on a very simple form: 



Mp)^{p) 



2i//(3i'-l) 



P 



1.539 



(8) 



As shown in Fig ^, this expression works remarkably well 
for larger densities into the semi-dilute regime. Devia- 
tions do occur for the dilute regime where Eq. (||) over- 
estimates the surface tension by a factor 1.5 for p 0, 
as demonstrated in the inset of Fig. |^. 

In a recent publication Maassen, Eisenriegler, and 
Bringer [ p9[ have used the renormalized tree approxima- 
tion to derive a surface tension which compares well with 
our results, as shown in Fig. |3| A similar asymptotic RG 
e-expansion compares slightly less well. The difference 
between the two approximations gives an estimate of the 
error in the RG approach. It should be kept in mind that 
our simulation approach also incurs small errors through 
the use of the fitted form of r(p), and because we use 
polymers of a finite length. 

Fuchs and Schweitzer recently applied the poly- 
mer reference interaction site model (PRISM) approach 
to polymer-colloid mixtures. In the limit of low colloid 
density, a number of analytic results can be derived for 
the insertion free energy Fi , from which the surface ten- 
sion can be extracted by using Eq. (|) and Eq. (||): 



„PRISM 



(p) = 1.279pRgil 



1.06—). 

p* 



(9) 



Here we have used the PRISM results arising from local 
packing information (see |^ for details) . As can be seen 
in Fig. H PRISM shows the correct qualitative but not 
the correct quantitative behavior. This is not surpris- 
ing since these results are based on a simplified PRISM 
model which has the advantage of being analytically solv- 
able, but the disadvantage of exhibiting the wrong scaling 
behavior. It would be interesting to see how well PRISM 
compares for the numerically more complex case where 
correct scaling behavior is included from the outset. 



III. DENSITY PROFILES AND SURFACE 
TENSION AROUND A HARD SPHERE 

Having described the surface tension for a polymer so- 
lution near a single hard wall, we next investigate the 
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related problem of a polymer solution near a single hard 
sphere (HS) of radius Rc- As discussed in the Introduc- 
tion, adding a single HS to a polymer solution reduces the 
number of configurations available to the polymers, and 
results in a finite insertion free energy or chemical po- 
tential described by Eq. (|l|). Besides the configurations 
directly excluded by the sphere of volume ^nR^, there 
are also configurations excluded near the surface of the 
sphere, an effect which manifests itself in an entropically 
driven depletion layer, just as was found for the case of a 
hard wall. However, the curvature of the sphere leads to 
corrections to the planar surface tension, as described by 
Eq. (^), i.e. the surface tension jsip) depends not only on 
the polymer density p, but also on R^ through the ratio 

q - Rg/Rc. 



A. Ideal polymers 

The free energy cost of inserting a single HS into a 
bath of ideal polymers is know Pq,E3 : 



(10) 



By combining this result with Eq. (|^) , it follows that the 
ideal polymer surface tension takes the form: 



iTip) = ll^ip) + pRgq 



The curvature corrections defined in Eq. ^ take on a 
particularly simple form here, since K'f{p) — pRg, and 
Ki ~ for i > 2. Note that this expression is not simply 
an expansion in q = Rg/ Rc] it is valid for all size ratios. 

In 1958 Asakura and Ooswawa ||l^ introduced a model 
where the ideal polymers are approximated as inter- 
penetrable spheres of radius Rao- This corresponds 
to approximating the true depletion layer by a step- 
function. The free-energy of insertion of a single hard 
sphere into a bath of AO particles can be easily calcu- 
lated to be 
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FC"" ^VAO—{l + qAo) 
Tao 



(12) 



where tjao = ^^P^'ao analogous to p/ p* , and we have 
defined the size ratio qAo = Rao/Rc- The surface ten- 
sion is therefore given by 

it'^ip) = it^'ip) + pRaoQao + ^ {qAof (13) 



where 7^*^ {p) — pRao ■ In this case the curvature correc- 
tions have a very simple geometrical origin: The volume 
of a spherical shell of width Rao with an inner radius of 
Rc has a larger volume than that of a flat layer of width 
Rao and area AttR^. In part this is a matter of defini- 
tion. For hard particles one can also find instances in the 
literature where Rc + Rao is taken as the radius of the 



Gibbs dividing surface. The AO model surface tension 
vanishes if it is defined in this way. 

If the prescription Rao = '^Rg/V^ is used to set the 
free parameter in the AO model, then the surface tensions 
for the planar wall is the same as that of ideal polymers. 
However, this prescription no longer holds for spheres im- 
mersed in a polymer solution since the curvature correc- 
tions to the surface tension for ideal polymers are not the 
same as those of the AO model. Physically this difference 
arises because the AO model assumes a fixed depletion 
layer width Rao while the (ideal) polymers can deform 
around a sphere, which leads to an effectively smaller 
depletion layer. This effect becomes progressively more 
pronounced with decreasing sphere size Rc/Rg Q. An 
effective AO parameter which takes this deformation ef- 
fect into account can be derived by equating the two sur- 
face tensions, Eq. (O) and Eq. (O) 



-"-AO 
Ra 




(14) 



Since the pressures in the two systems are the same, i.e. 
n = p, this is equivalent to equating the two insertion 
free-energies Fi of Eq. ( p^ ) and Eq. (|l2|), as done in ref- 
erences [^^. For q = Rg/Rc this expression re- 
duces to R'^/q I Rg ~ 2/ the known result for a single 
wall. For large q on the other hand hand, the effective 



(11) AO radius scales as R'^Q/Rg ~ q For fixed Rg, 

the effective radius R'^/q decreases monotonically with 
decreasing sphere size Rc, as shown in Fig. 0. 



B. Interacting polymers 

1. Low density limit for interacting polymers 

For interacting polymers, the p —^ limit of the 
curvature corrections to the surface tension have been 
calculated to first order in an e-expansion by Hanke, 
Eisenriegler, and Dietrich For large spheres (small 
Rg/Rc), they find: 



lim 



Isip) 



1 + 0.849? 
0.0375g2 + 0{q^) 



(15) 



In this low density and small q regime, the curvature 
corrections for interacting polymers are quite similar to 
those found for non-interacting polymers. Compare, for 
example, the first relative curvature correction coeffi- 
cient, which is 0.849 for interacting polymers, and 0.886 
for ideal polymers. In the opposite (small sphere) q 00 
limit the differences are more pronounced: 7s'^/7 



id 



AO UAO 



ilAo)^ while for interacting polymers RG 



and scaling theory approaches predict that 7s /7« 

(g)l/--l « g0.70 



1 [|l]j^ 



4 



2. Interacting polymers at finite densities 

We have calculated the density profiles h{r) = p{r)/p— 
1 for polymers around spheres of radius Rc — lA9Rg, 
Rc = 0.95Rg and Rc = 0.59i?c from simulations of 
L — 500 SAW polymers. These are shown in Fig H in 
the CM representation. Just as was found for the case of 
a planar wall, the depletion layers shrink with increasing 
bulk density p. Because the polymers can deform around 
the colloid, the density profiles in the CM representation 
can penetrate into the HS region, an effect which becomes 
more pronounced for smaller colloids (larger q). (For an 
interesting proposal that describes the monomer density 
around a spherical particle we refer to reference fi^ .) 

The relative adsorption around a sphere is defined as: 



AirRltsip) 



p dp 



Anr^ h{r)dr + -nR^ 
3 

(16) 



Here h{r) is defined from the center of the sphere, and 
^sip) has the dimension of a length. The volume of a 
the single HS was subtracted off so that the adsorption 
only describes the effects of the depletion layer around 
a sphere. For low density the relative adsorption of a 
sphere is larger than that of a planar wall by a curva- 
ture correction factor term similar to those described in 
Eq. ( p^ ) . As the density increases the relative adsorption 
decreases and tends asymptotically to the same value as 
for a planar wall. This can be understood from the sim- 
ple "blob" picture in the semi-dilute regime: Since 
the ratio of the blob-size to the sphere ^ (p) /Rc decreases 
with increasing density, the curvature corrections to the 
relative adsorption are also expected to become relatively 
less important with increasing density. 

The surface tension 7s (p) can now be calculated from 
Eq. (^) using the adsorption defined in Eq. ([T6|). In Fig.||, 
we compare the surface tension for three different sphere 
sizes to 7to(p), the value for a planar wall. As expected 
from the results for low densities, (see e.g. Eq. (|lq)), for 
a given density p, the surface tension increases with de- 
creasing Rc- The ratio ")s[p)hw{p)^ shown in Fig. de- 
creases with increasing p/ p* . Again, the rate of change 
is largest in the dilute regime; for increasing p/ p* the 
two terms appear to approach each other asymptoti- 
cally. Just as we argued for the adsorptions, the "blob" 
picture in the semi-dilute regime implies that the cur- 
vature corrections should decrease with increasing den- 
sity, which is what we observe. This also implies that 
ls{p) ~ lw{p) ^ pi-539 ^]-^g semi-dilute regime. Of 
course the smaller the HS, the higher the density one 
needs for the curvature corrections to become negligi- 
ble. This picture is confirmed by recent scaling and RG 
arguments |29[ , which show that the first curvature cor- 



rection coefficient ki ~ p 



„0.770 



implying that 



with increasing density, the contribution of the curvature 
corrections defined in Eq. (0) becomes relatively smaller. 



so that "fs approaches 7u, . In the inset of Fig. |^ we com- 
pare our results to the RG calculations, valid to lowest 
order in q, i.e. 7s (p) — ^w{p) + Ki{p)q. Although only 
the results for the ratio q — 1.05 are shown, they are 
similar to those at the other two size-ratios, which also 
show an overestimate by the RG. The difference may be 
due in part to higher order Ki{p) terms which have not 
yet been calculated by RG. To confirm this picture fur- 
ther simulations are needed since: (a) Our simulations of 
the adsorption are only for p/ p* < 2.32, and we extrapo- 
lated to higher densities using a fit form which scales as 
-fs(p) ^(p) p'0.770 ^-^ j^jgj^ densities, (b) We only 
examined 3 different sphere sizes so that it is difficult 
to directly extract Ki{p), and for that matter the higher 
order Ki{p). 

Finally, we reemphasize how much the density de- 
pendence of the surface tension of the interacting poly- 
mers differs from that of ideal polymers or the related 
Asakura-Oosawa model, where the ratio of the wall to 
the sphere surface tensions is independent of density, and 
close to that of interacting polymers in the low-density 
limit (Compare Eqs. ([l^), (|l|), and (|l|)). 

3. Comparison with a hard sphere model 

One might inquire what would happen if the polymers 
were modeled as HS instead. By using the very accurate 
Rosenfeld fundamental measure density functional |4^ ] 
technique, an explicit form for the surface tension of a 
HS fluid with radius Rp around a single HS (radius Rc) 
has been calculated H] 



Is (rip) 



= 1 



2(1 



(2- 



2{1 - rjp)Hnil ~ rjp) ^Rp 
3r]p{2 + r]p) Rc 



(17) 



where we have defined the packing fraction rjp — ^TrppR^, 
for a number density pp. The planar wall-surface tension 
is given by 



IwiVp) 



3r]p{2 + r]p) 
87ri?2(l - rjpy 



(18) 



We note that this result for the surface tension of a HS 
fiuid around a sphere was also derived independently by 
scaled particle theory |45|. Eq ( [l7| ) can be generalized to 
the non-additive HS model, for which the cross diameter 
Rep 7^ ^{Rc + Rp)^ so that one can smoothly interpo- 
late between the fully repulsive HS model and the fully 
non-additive AO model 

To lowest order in density, the surface tension of a HS 
system near a planar wall is 7^'^ ~ RpP, i-e. the same as 
that of the AO model, which is close to that of interacting 
polymers in the same limit where ~ 1.074i?gp. How- 
ever, the terms of order p^ are already significantly larger 
in the HS case. Therefore, as illustrated in the inset of 
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Fig H, the HS model gives a large relative overestimate 
of the surface tension well before reaching the packing 
fraction at which the system freezes, (here we took units 
where Rp = Rg so that 77^ is equivalent to pj p*\ 

For a fixed size-ratio Rp/Rc the curvature corrections 
for a HS system vary with density as: 



AO 



where the ratio for the AO model comes from Eq. (|l|) 
with Rao = Rp- As illustrated for a 1:1 size ratio in the 
inset of Fig. ^, for small -qp this ratio is indeed almost lin- 
ear. The change with density is more pronounced than 
that found for a polymer-colloid system with a similar 
size-ratio, suggesting (not surprisingly) that a full HS 
system is not such a good model of interacting polymers, 
even at relatively low densities. Making the spheres non- 
additive does not fundamentally alter this picture - the 
behavior of polymers falls into the class of "mean field 
fluids" i.e. they do not behave like hard-core par- 

ticles. 



4- Direct calculation of Fi by the Widom insertion method 

We also performed direct computer simulations of the 
free energy Fi by measuring the insertion probability 
of a single sphere in a bath of polymers at fixed den- 
sity p/p* = 1.16. This is closely related to the so-called 
Widom insertion technique to find the chemical potential 
psf . Fig. H shows that Fi grows with increasing sphere 
size as expected. The same is true for the contribution 
due to the depletion layer, i.e. the contribution propor- 
tional to AttJsRc ^1- (0)- However, the relative impor- 
tance of this surface tension term increases with decreas- 
ing sphere size, and becomes the dominant contribution 
as Rc/Rg — > 0. The values up to Rc/Rg = 0.59 were 
calculated by the insertion probability method, while 
those with larger Rc/Rg were taken from the adsorption 
method, i.e. from the density profiles, as was done for 
example in Fig. ^. For Rc/Rg — 0.59 we used both meth- 
ods and find very similar results, suggesting that the two 
approaches are mutually consistent. We also compare to 
results for ideal polymers (Eq. and for PRISM [O. 



5. Limit of small colloids 

In the limit of small Rc/Rg, scaling arguments and 
RG theories predict that the free energy to insert a sin- 

t particle in a bulk polymer solution takes the form 



Fi 



AgRc 



(20) 



Where Ag is a universal numerical pre-factor that can 



Eq. ^ reduces to Fi « 18.4pi?i-30i?i-70. This expres- 
sion is directly compared to our simulations in Fig. g|. By 
comparing to Eq. |l| we can extract the surface tension 
from the insertion free energies. This was done for the F 
at /9//0* = 1.16 shown Fig. @, and also for L — 2000 poly- 
mers at p/p* ~ 0.94. Using the longer polymers allows 
effectively smaller colloidal RcS to be used in our lattice 
simulations. The surface tensions are depicted in Fig. ||. 
At small q our computer simulation results correspond 
reasonably well with the asymptotic RG results. We ex- 
pect there to be small errors due to the discreteness of 
out lattice simulations, similar to those depicted for ideal 
polymers in Fig. ^. These discretization errors become 
more important as the spheres become relatively smaller. 
We have made some small corrections j49| to take this 
into account, but a more systematic study, possibly with 
longer polymers, would be necessary to completely test 
the RG results. 

When the colloids are much smaller than the polymers, 
one expects that they only probe the local monomer 
density, and not the overall number density of poly- 
mer coils. In fact, Eq. (20) implies just that since 
Fi cx pR^^Rl-^" - pLRl-^" = cRl-^°. The reason Fi 
scales linearly with the monomer density c is that by 
definition this is very small (c <C 1) in the dilute and 
semi-dilute regime. The small colloidal particles probe 
what is effectively an ideal gas of monomers. 

For ideal polymers Fi cx cRc in the limit of small 
i?c, which implies that for a given Rg and p, and for 
a small enough Rc it is easier to insert a hard-sphere into 
a bath of interacting polymers than it is to insert it into 
a bath of non- interacting polymers. At first sight this 
may seem surprising, but the reason is as follows: In- 
side an interacting polymer, the monomer concentration 
scales as c ~ i? 
c - 

mer and make it less dense; for a given Rg, the monomer 
density c is larger for ideal polymers than for interaction 
polymers, and since the small colloids only probe the lo- 
cal monomer density it is easier to insert the sphere into 
an interacting system than into a non-interacting system 
at the same Rg. This effect is illustrated in Fig. for 



^-1.30 for ideal polymers it scales as 

^. In other words, the interactions swell a poly- 



P/P* 



1.16, where the crossover is at about Rc ~ O.bRo 



be calculated from an RG technique 



For d = 3 



(This limit should not be confused with a comparison at 
fixed L). Note that the PRISM results also overestimate 
Fi at small Rc- This is in part because the simplified 
PRISM model we compare to also includes ideal polymer 
statistics, resulting in an overestimate of the monomer 
density compared to a true interacting system, an effect 
already pointed out in ref. [^ . 

For large spheres, on the other hand, where Fi w 
^nIl{p)Rc which scales as Fi ^ p^-^^R^ in the semi- 
dilute regime, the spheres do directly probe the number 
density of polymer coils, and the insertion free energy for 
interacting polymers is always higher than that of ideal 
polymers at the same Rg and p. Note how differently 
the large and small Rc/Rg hmits of Fi scale both with p 
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and with i?c- Significant differences can be also found for 
the scaling of the surface tensions since for large Rc/Rg, 

Is - -1M^ + 0{RJR,)) ^ pl-539(l + O(p-0-770i?^/i?^), 

while for small Rc/Rg, the RG expressions imply that 
Is - cR-° '. 



C. Surface tension for polymers as soft colloids 

We have recently modeled polymers as single "soft col- 
loids" interacting with a pair potential between their 
CM Hfll. These pair potentials were derived by a 
liquid state theory based inversion procedure such that 
the soft colloids have exactly the same radial distribu- 
tion function g(r) as those generated by a fully interact- 
ing polymer simulation. A similar inversion procedure 
was used to derive the potential between the soft-colloids 
and a planar wall or a HS. These wall-polymer or sphere- 
polymer potentials are such that they exactly reproduce 
the one-body density profiles p{r). 

Since our effective polymer-polymer potentials pro- 
vide a very accurate representation of the pressure n(p) 
p9| , pl| , while the polymer- wall or polymer-sphere inter- 
actions are constrained to reproduce the correct den- 
sity profiles, and therefore the correct adsorption T{p), 
Eq. (||) implies that our soft-colloid approach has the 
correct surface tensions automatically built in. Similarly 
Eq. implies that this approach correctly reproduces 
Fi for a sphere immersed in a polymer solution. 



IV. CONCLUSIONS 

In summary then, we have used computer simulations 
of SAW polymers on a cubic lattice to calculate the den- 
sity profiles for non-adsorbing polymers near a planar 
wall, and near HS's. From this we were able to calculate 
and fit the relative adsorption r{p). Together with the 
equation of state, which is well understood for polymer 
solutions, this provides the needed ingredients to calcu- 
late the surface tensions through Eq. (g). 

The surface tension of interacting polymers near a 
planar wall was shown to differ significantly from that 
of ideal polymers, or other simple models such as the 
Asakura Oosawa penetrable-sphere model, or a pure HS 
fluid. Similarly, a recent PRISM calculation also 
shows large qualitative differences with our results, which 
could have been anticipated in view of its use of simpli- 
fied ideal polymer statistics. On the other hand, some 
recent RG results compare very well to our calcula- 
tions. In the semi-dilute regime, the surface tension sim- 
plifies to the form given in Eq. (H), which implies that 

Near a sphere with a radius of the same order or larger 
than Rg, the surface tension 7s (p) of the polymer solution 
can be written in an expansion in the size-ratio q. For de- 
creasing sphere size (increasing q), the ratio 7s(p)/7t«(p) 



increases. For a given q, however, 7s(p) approaches jwip) 
as the density increases. We attribute this to the decrease 
of the effective curvature corrections with increasing den- 
sity since the blob size scales as (,{p) ~ p-o.770 ^j^g 
semi-dilute regime. This is again consistent with some 
recent RG calculations of the correction coefficient ki (p) 
, although further simulations are needed to confirm 
the scaling and form of the coefficients Ki{p). 

For smaller colloids, it is advantageous to use a direct 
Widom insertion technique to calculate the free-energy 
-Fi. For very small colloids (large q), our simulations 
were consistent with asymptotic RG results which sug- 
gest that Fi oc cR^'^°. This insertion free-energy is dom- 
inated by the contribution of the depletion layer; its be- 
havior is qualitatively different from the behavior found 
at smaller q, and suggests that the expansion of Eq. (|^) 
breaks down for large q. 

Because our soft-colloid approach was derived to repro- 
duce the correct one-particle density profiles near hard 
walls or hard-spheres, it will automatically reproduce the 
correct adsorptions, and therefore also the correct surface 
tensions and insertion free-energies -Fi . 

The walls and spheres in this study are purely repul- 
sive. Adding a wall-polymer attraction should decrease 
the amount of depletion, and therefore also lower the sur- 
face tensions. More subtle effects could be expected if in 
addition the solvent quality is decreased. New effects are 
also expected for binary mixtures of polymers with selec- 
tive adsorption of one of the species. These systems will 
be the subject of future investigations. 

The next step is to move from the one-sphere or one- 
wall problem to the case of a two-sphere or a two-wall 
system, and calculate the effective interactions between 
the two particles. This is the subject of a forthcoming 
paper [pO| - 
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FIG. 1. (a) The wall-polymer CM profile 

h{z) = p{z)/p — 1 for L — 500 SAW polymers at different 
bulk concentrations, (b) The wall-monomer profile hm{z) for 
the same bulk concentrations. Both representations result, by 
definition, in the same relative adsorptions. The full lines are 
a fit to the simple form hm{z) = tanh^(2;/r(p)) — 1 




FIG. 3. Polymer- wall surface tension "fw{p) divided by 
7^(p). The full lines with square symbols are for interacting 
polymers and were calculated with Eq. (^), while the dot- 
ted line with circle symbols denotes the simpler expression of 
Eq. (^) which is only valid in the semi-dilute scaling regime, 
where 7^, ~ p^'^^^. Also shown are two recent RG results 
[ p9| : the dashed line denotes the renormalized tree expan- 
sion, while the dot-dashed line denotes the asymptotic limit 
for an e-expansion. The solid line is from Eq. (H) , a result de- 
rived from a recent PRISM calculation |Q. The inset shows 
the ratio of the full and simplified expressions for 7™(p). They 
coincide for higher densities but in the low density limit, the 
semi-dilute scaling expression overestimates the true surface 
tension by a factor 1.5. 
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FIG. 2. Relative adsorption f(p), in units of Rg, ver- 
sus density. Circles denote direct computer simulations of 
L = 500 SAW polymers near a single hard wall, and the line 
denotes the simple fit with the correct scaling behavior, given 
by Eq. (^). Also shown is the density dependence of the ra- 
dius of gyration. In the semi-dilute regime T « — ^ ~ p°''^™, 
while 7?g(p)/i?g(p = 0)~p« "5 
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FIG. 4. The effective AO radius R^/^, given by Eq. (0), 
which would result in the same surface tension for an AO fluid 
around a sphere of radius Rc as found for ideal polymers of 
size Tig. For infinite sphere size (i.e. a wall) Rao = {'^/^/^)Rg- 
As the relative sphere size Rc/Rg decreases this effective pa- 
rameter decreases due to the deformation of the polymers 
around a sphere. The symbols denote direct simulations 
L = 200 ideal polymers on a lattice, taken from reference 
[Q. The small differences are due to the discrete nature of 
the lattice used in the simulations. 
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FIG. 6. Surface tension for a planar wall, and for spheres 
with q = 0.67, q = 1.05 and q — 1.68 as a function of density. 
We also include the planar surface tension of a HS fluid, with 
Rp = Rg such that rjp = p/p* . Inset: Blowup of the graph for 
low densities. The planar wall surface tension for interacting 
polymers, ideal polymers, and the HS system are compared. 




FIG. 5. The polymer CM density profile h{r) around a 
sphere for the ratios q = Rg/Rc — 0.67,1.05,1.68 (graphs 
from top to bottom). For each sphere size the curves are 
for p/p* = 0.037,0.30,0.59,1.17,2.33 (solid, dotted, dashed, 
long-dashed, and dot-dashed lines respectively). The deple- 
tion layer narrows with increasing density, just as was found 
for as single wall (compare with Fig. The small vertical 
lines denote the position of the radius of the colloid. The 
polymers can wrap more easily around the smaller colloids, 
which explains why the CM profile penetrates further into the 
colloid for smaller Rc/Rg. 




FIG. 7. Ratio of surface tension of a sphere to the sur- 
face tension of a wall for spheres with q — 0.67, q — 1.05 and 
q — 1.68. Inset: Comparison of an RG calculation valid 
to lowest order in q, and our direct calculation for q — 1.05. 
We also compare the ratio of the surface tension of a HS fluid 
around a single inserted sphere to the planar HS surface ten- 
sion. The size-ratio is 1 : 1, and rip = p/p*. The value at 
77p = is equal to that of the AO model, given by Eq. (|l|). 
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FIG. 8. Insertion free energy Fi for spheres of various radii 
Rc, in a polymer bath at p/p* = 1.16. For smaller Rc/Rg a 
direct insertion method was used, while for larger Rc/Rg the 
adsorption method was used. We also compare 47r ii^ 7s (p), 
the contribution to Fi due to the creation of a depletion layer. 
For small Rc/Rg this term is the dominant contribution to the 
insertion free energy Fi. Comparison is also made to an ex- 
pression from RG theory, Eq. (po|), valid in the small Rc/Rg 
limit with results from the PRISM approach |^ and 

with Fi for ideal polymers, taken from Eq. (p^). Note that 
for this density, the ideal and interacting results for Fi cross 
each other at Rc/Rg ~ 0.5, below which it is easier to insert 
a spherical colloid into an interacting polymer solution than 
into a non-interacting one. 
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FIG. 9. Surface tension 7^ for spheres of different radius 
Rc, in a polymer bath at p/p* = 1.16, (L = 500 SAW simula- 
tions) and for p/p* = 0.94 (L = 2000 SAW simulations). The 
insertion and the adsorption methods agree to within the ex- 
pected statistical errors of our approach for Rc — 0.59Rg. We 
also compare to an expression from RG theory, Eq. ( po| ) , valid 
in the small Rc/Rg limit |^^. The arrows on the right depict 
the values of the planar surface tensions, valid as Rc/ Rg ^ 00. 
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